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■ Abstract 

■ A population has two types of individuals, each occupying an island. 

One of those, where individuals of type 1 live, offers a variable environ- 
ment. Type 2 individuals dwell on the other island, in a constant en- 
vironment. Only one-way migration (1 — >■ 2) is possible. We study the 
asymptotics of the survival probability in critical and subcritical cases. 



>' 

\0 '. 1 Introduction 

I We study a two-type branching process in random environments. It can be 

■ viewed as a stochastic model for the sizes of a geographically structured popu- 

0^ I lation occupying two islands. Time is assumed discrete, so that one unit of time 

i represents a generation of individuals, some living on island 1 and others on 

^ I island 2. Those on island 1 give birth under influence of a randomly changing 

environment. They may migrate to island 2 immediately after birth, with a 
probability again depending upon the current environmental state. Individu- 
. als on island 2 do not migrate and their reproduction law is not influenced by 

^ ■ any changing environment. Our main concern is the survival probability of the 

whole population. 

An alternative interpretation of the model under study might be a population 
(type 1) subject to a changing environment, say in the form of a predator 
population of stationary but variable size. Its individuals may mutate into a 
second type, no longer exposed to the environmental variation (the predators, 
e. g.). 

The model framework is that furnished by Bienayme-Galton- Watson (BGW) 
processes with individuals living one unit of time and replaced by random num- 
bers of offspring which are conditionally independent given the current state of 
the environment. We refer to such individuals as particles in order to emphasize 
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the simplicity of their hves. Particles of type 1 and 2 are distinduished according 
to the island number they are occupying at the moment of observation. Our 
main assumptions are: 

• particles of type 1 form a critical or subcritical branching process in a 
random environment, 

• particles of type 2 form a critical branching process which is independent 
of the environment. 

In Section [2] we recall known facts for constant environments. They will 
then be compared to the results of this paper on random environments. In Sec- 
tion [37T] we describe IID environments (Independent and Identically Distributed 
environmental states), and then in Section [321 Markovian environments. 

Our first result on the asymptotics of the survival probability, Theorem [T] 
is stated and proved in Section It deals with the critical case with an IID 
environment. Section [5] contains our result for the subcritical case and an IID 
environment. The last two Sections, [6] and [3 expand our results for the IID 
case to the Markovian setting. 

Use of notation. In asymptotic formulae constants denoted by the 
same letter c are always assumed to be fixed, that is independent 
of the parameter that tends to infinity (or zero) . 

2 Two-type decomposable branching processes 

Consider a two- type BGW-process initiated at time zero by a single individual of 
type 1 . In this paper we focus on the decomposable case where type 1 particles 
may produce particles of types 1 and 2 while the type 2 particles can give birth 
only to type 2 particles. Put 

• Xn = the number of particles of type 1 present at time n, in particular, 

Xo = 1, 

m Yn — the number of type 2 daughters produced by Xn particles of type 1, 

• T = the first time n when Xn = 0, 

• Z„ = the number of particles of type 2 present at time n, as a rule, we 
assume that Zq = and, as a result, Zi ~ Yq, 

• Sn = X]fc=o -^f' ^^^^ gives the total number ever of particles of type 
1. 

• Wn = J2k=o so that Wt gives the total number of type 2 daughters 
produced by all St particles of type 1. 
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The aim of this section is to summarize what is already known about such 
branching processes in the case of a constant environment. This will pave our 
way in terms of notation and basic manipulation with generating functions 
towards branching processes in IID random and then Markovian environments. 

If the environment is constant from generation to generation, two-type de- 
composable BGW-processes are fully described by a pair of probability gener- 
ating functions 

h{s)=E[s^], 

where and ^2 represent the numbers of daughters of type 1 and 2 of a mother 
of type 1, while 77 stands for the number of daughters (necessarily of type 2) of 
a mother of type 2. Let 

„ - IT r/r 1 - I 

fll —tj [4iJ — — \s^=S2=l ) 

1^2= - J^jJ = Ui=S2=l ) 

t>l-tj [42J - Ul=S2=l I 

C^2 = E [^2(^2 - IjJ = Ul=S2 = l > 

mi=E[r,] = h'il), 

m2 = E [77(7? - 1)] = h"{l), 

be the first two moments of the reproduction laws. Concerning the second type 
of particles we will assume that 

TOi = 1, 777-2 G (0,00), (1) 

implying that the probability of extinction 

Q„ = P(Z„ = 0|Xo = 0,Zo = l) 

(of a single-type BGW-proccss evolving in constant environment with the prob- 
ability generating function h{s)) satisfies 

2 

i-Qn , 71 -I- 00. (2) 

771271 

It follows that 

an ■■= -l0g/(l,(9„) ~ 1 - f{l,Qn) ~ , 71 00. (3) 

771271 

We will be interested in two kinds of reproduction regimes for particles of type 1, 
critical and subcritical. In the constant environment setting with fi2 € (0,oo), 
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the critical case corresponds to /xi = 1 and the subcritical case is given by 
/ii S (0, 1). In the critical case with a constant environment we have 



P (Xn > 0) = P (T > ?i) , n^oo, 

H2n 



and according to Theorem 1 of 



P (X„ + Z„ > 0) - P (Z„ > 0) 



2v^ 



, n — >■ oo. 



Next we outline a proof of ([S|) based on the representation 



(4) 



(5) 



P (Z„ > 0) = E 



E 



fc=0 



E 



1_ ■Q/^.(1^Q^_^) 



fe=0 



2 _ g- 5i;fc=o -^fca„_fc 



(6) 



preparing for the proof in the random environment case, to be given in Section 
HI Thanks to (g]) and 

P (Z„ > 0) < P (X„ + Z„ > 0) < P {Xn > 0) + P (Z„ > 0) (7) 

in order to prove ([5]) it is enough to verify that 

P iZn > 0) ^ ■ n -5> CXD. 

However, by the branching property the total progeny of a single-type branching 
process St is 1 plus independent daughter copies of St- In terms of the 
Laplace transform 

0(A) =e-V(0(A),l), 
where 0(A) ~ Ee^'^'^^. As A — > Taylor expansion yields 

1 - 0(A) = 1 - e-^ + e-Vi(l - 0(A)) - £-^^(1 - 0(A))2(1 + o(l)). (8) 

For /^i = 1, this implies - after some calulation - that 

1 - 0(A) ^ y^2}:/jr2, A ^ 0. 

Thus, due to © 



E [1 



and it remains to verify, see ©, that 



E 



1 - e 



E [l - e-^^"^"] , 
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This holds, indeed, since by (g]) and for any fixed e > the probabihty P(T > 
ne) is much smaller than the target value of order c/y/n. (In [19] and [20] infinite 
second moments in decomposable two- type critical processe were allowed.) 
On the other hand, in the subcritical case ([H]) implies that 

1 - 0(A) -A/(l- Ail), A^O, 

so that by © 

7712(1 ~fii)n 

In view of P (X„ > 0) ^ c/i" we conclude that in the subcritical case 

26*1 

P(X„ + Z„ >0)-P(Z„ >0) — , n^c». (9) 

m2(l - ^i)n 

See [15j for a comprehensive study of subcritical decomposable branching pro- 
cesses in a constant environment. 

3 Branching processes in a random environment 

A randomly changing environment for BGW-processes is modeled by a random 
sequence of probability generating functions for the offspring distributions of 
consecutive generations. Throughout this paper we assume that the offspring 
distribution for type 2 particles is the same across the different states of the 
environment and characterized by the same generating function h{s). This 
restriction greatly simplifies analysis still allowing new interesting asymptotic 
regimes. 

We consider two types of stationarily changing environments: IID and Marko- 
vian. 

3.1 IID environment 

We start our description of the IID environment case by a simple illustration 
based on just two alternative bivariate generating functions f^^\si,S2) and 
/^^^(si, S2) with mean offspring numbers (/i^^-* , 0^^'') and {fj,'f\9'^'') respectively. 
We assume that at each time 77 the environment is say "good" with probability 
TTi, so that the type 1 particles reproduce independently according /^^■'(si, S2), 
and with probability 112 — 1 — tti the environment is "bad" and particles of 
type 1 reproduce according to the /'^■'(si, S2) Ie^w. In other words, the generat- 
ing function f {81,82) should be treated as a random function taking the form 
/^^••(si, 82) or /^^'(si, 82) with probabilities tti and 7r2. In particular, the vector 
of the mean offspring numbers (/ii,6'i) takes values (/i^^',0^^') and {fj,'^\9^^) 
with probabilities tti and 7r2. 

More generally, our two-type branching process in an IID random environ- 
ment is characterized (besides the fixed reproduction law ft,(s) for the type 2 
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particles) by a sequence of generating functions {/n(si, S2)}$?Lo independently 
drawn from a certain distribution over probability generating functions so that 

/«(S1,S2) = /(si,S2). (10) 

In this setting the respective conditional moments /ii, /i2, di, and 62 should 
be treated as random variables. An important role is played by the random 
variable C = log fii representing the step size of the so-called associated random 
walk [3] formed by the partial sums Co + ■ ■ ■ + Cn-i with Ci = C- 

Use of notation. Characteristics of the reproduction law in gen- 
eration k are denoted by adding an extra lower index k to the 
generic notation, like in (ITUl) . Another example: {Oi^ki(^2.k) in (11) 
ipSl) represent the conditional moments {61,62) for the offspring 
numbers of the type 1 particles living at time k. 

3.2 Markovian environment 

One way to relax the IID assumption on the environment is to allow for Marko- 
vian dependence among its consecutive states. We implement this idea by 
modeling changes in terms of an irreducible positively recurrent Markov chain 
{cnl^o with countably many states {1, 2, . . .}. Assuming a stationary initial 
distribution, we will associate with each state i of this chain a probability gener- 
ating function (si, S2), so that the changing environment for the branching 
process is governed by the sequence of identically reproduction laws 

/n(si, S2) = f'-''"\si,S2), n = 0, 1, . . . 

with Markovian dependence. Due to the stationarity we can again write (jlOp 
and use the same notation for the marginal moments of the reproduction laws 
as in the IID case. 

To build a bridge to the IID environment case we use an embedding through 
a sequence of regeneration moments {t^I^q defined as 

To := 0, Tfc+i min{n > : e„ = Cq}. (12) 

The times Tk+i — Tk between consecutive regenerations are independent and all 
distributed as r ti. The embedded process {X*, Z*) defined as 

iX:,Z:) = iX^^,ZrJ, n = 0,l,... 

is a decomposable branching process in an IID environment with two types of 
particles 1* and 2* and conditional reproduction generating functions 

r (si, S2) /(^«) (/(^^) ( . . . {f^'-'\s,,S2), h{s2)) ...), K-1{S2)) , (13) 
h*{s) := h{h{. . . h{s) ...)) = hr{s), (14) 

where /ifc(s) stands for the fc-fold iteration of h{s). 
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The key difference from the IID case is that the reproduction law for the 
2*-type particles being dependent on the random environment. However, this 
dependence is of specific nature which we are able to manage using the law 
of large numbers for the renewal process. Notice that on its own the 2*-type 
particles form a so-called degenerate critical branching process in an IID ran- 
dom environment [3]: its conditional offspring mean is deterministic m| — 1. 
Meanwhile, the conditional variance is random mi, = Tm2. 

Taking the first and second order derivatives of ([13]), we can express the the 
moments of the star-reproduction law in terms of the moments of the consecutive 
reproduction laws with Markovian dependence. In what follows we use (jlll) 
again, while keeping in mind that the sequence {fJ-i^k, l^2,k,9i_k,02.k)kZQ no'^ 
consists of dependent random vectors. It can be shown that 

r— 1 r — 1 r — 1 

iJ'*i=Y[fj'i,k, f^*2 ^ fj-*!^— n '^i^*' 

fc=0 fc=0 ^^''^ i=k+l 

T-1 fc-1 

^1 = X] ^i-'" n ^1'*' 

k=0 i=0 



Furthermore, setting 



we can write 

T-1 k-1 



n j-1 
Ak,n = ^ Ol.j Y\_ A^l.* 



j=k 



k=0 i=0 

r-2 k-1 

(/U2,fe-4^+l,r-l + '^^J-lM^lM^k+l.T-l + o-^ (r - 1 - fc) 9i^k) Y\_ A^i/ 
fc=0 i=0 



Considering 



T-1 



fc=0 

observe that due to stationarity of the Markovian environment we can use a 
version of the Wald identity 

Er]=E[r]E[C] (15) 

provided 

oo 

^E[|a|l{T>fe+i}] <oo. (16) 

k=0 
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Equation follows from 

oo k—1 oo 

E [CI = E E E [c«i{r=fe}] = E E [c(e„)i{.>„}] 

cxD oo oo 

= E E C(*)P(en = r > n) = ^ C(*)^^E [r] = E [r] E [C] , 

n— 2—1 2—1 

where tt^ are the stationary probabilities and the second last equality is justified 
as follows. According to Theorem 6.5.2 from [TO for any state j 

oo 

A'j(*) = E^(^" = i,T> n\eo = j) 

ri=0 

defines a stationary measure which is necessarily of the form = CjTTi, where 

oo oo oo oo oo oo 

E '^o^i = E E = E E E P^"^" = i, r > n|eo = j) 

J — 1 J — 1 i— 1 j/ — 1 n— i— 1 

oo 

= ^^,E[T|eo=j] = E[r]. 

Developing the example of two environmental states from Section 13.11 let us 
consider a Markov chain {e„}5iJ°^Q with transition probabilities 

/ 1 - dTT2 d-K2 ^ d min ( ^ M 

\^ dTTl 1 - dlTl y ' V TTl ' 7r2 / 

and a stationary distribution (771,712). (Notice that d = 1 corresponds to the 
IID case.) Assuming stationarity we get for the regeneration time 

P (t = 1) = 7ri(l - d7r2) + 7r2(l - din) = 1 - 27ri7r2d, 

P (r A:) = 7ri(i7r2(l - d7ri)''"^a7ri + 7r2d7ri(l - dTT2)''^^ dTr2 

= d7ri7r2(d7ri(l - dni)''-^ + d7r2(l ~ ^772)''"^) , fc > 2, 

implying 

E[r-l] = l, E[r(r-l)] = -^-i. 

dniTr2 d 

If (61, 62) are the two possible values for (, we can write E [(] = vri^i +7r262 and 

E [CI = E [E [C* |r] ; r = 1] + E [E [C*|r] ; r > 2] 
= 7ri(l — dn2)bi + 7r2(l - ^711)62 

+ TTldTT2 (^h + + 7r2d^l (^62 + = 2E [C] , 

confirming equality p^ . 
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4 Critical processes in IID environment 

The single type critical branching process with an IID environment displays an 
asymptotic behavior that is in stark contrast with the constant environment 
formula According to Theorem 1 in [T3], if 

E[C]-0, Var[C] G (0,oo), (17) 

E [^J.2^h^ (1 + max(0,log^i))] < cx), (18) 

then 

P(X„ > 0) = P(r > n) - " ^ oo. (19) 

(A much more general limit theorem is obtained in [5].) The following theorem 
shows that in the decomposable case the difference between the constant and 
random environments is even more striking. For constant environments the 
survival probability decays as c/-\/n, see (O, but in random environments the 
decay is like c/ logn. 

Theorem 1 Consider a critical decomposable branching process in an IID en- 
vironment satisfying ([I]), ( | j7| ), (ITS)) , and 

E [fi^^] < oo. (20) 

If for some positive a 

P {01 <l/x) = o {{logx)-^-") , x^oo, (21) 
P{ei>x)^o ((log x)-^-") , x^oo, (22) 
P(6i2 > a;6ii) =o((loga;)"^"") , a; ^ oo, (23) 

then there exists a constant Kq such that 

P(X„ + Z„>0)-P(Z„>0)--^, n^oo. (24) 

log n 

Before the proof, we make some comments on the conditions and statement 
of this theorem. 

Use of notation. From now on we will often use abbreviations 

Xa = (log a;)^+" and ria = (log n)^+". ^ ^ 

Conditions (|2T|) . ([22l) . and ([23]) are needed for the following properties to hold 
for any fixed e > 

P( min 01 k < x^"] = o( . X oo, (26) 

P f max 6*1 fc > I = o f ) , x oo, (27) 

\0<k<x^ ' J \\ogx J 

max {02,k/Oi,k) > xA =0 [ t— \ , x ^ oo. (28) 
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Each of them is proven via intermediate step Hke 



P ( min 6ik <x " ) < XaP {9i < x 

,0<k<Xa 



relying on the IID assumption for consecutive environmental states. The con- 
stant Kq in the statement (|24|) is the same as in the asymptotic formula from 

m 



P (St > x) 



Kq 

loffx 



, X — > oo. 



(29) 



This constant has a complicated nature and is not further explained here. It is 
necessary to mention that the representation ([^ has been proved in [1] under 
conditions P7|) , (|18p , and (I^U)) only for the case when the probability generating 
functions /n(s,l) are linear-fractional with probability 1. However, the latter 
restriction is easily removed using the results established later on for the general 
case in [13] and [3]. 

We now turn to the proof of Theorem [T] To derive (p4)) it suffices, in view 
of (Uni), to show that 



P {Zn > 0) 



Kq 
logn 



, n — > C!0, 



or 



liminf {logn • P(Z„ > 0)} > Kq, 
lim sup {logn • P(Z„ > 0)} < Kq. 



(30) 
(31) 



These inequalities will be established next using a counterpart of ([6]) 



P (Z„ > 0) = E 



E 



1-n 



■' n — k 



k=0 



E 



k=0 



1 - CXp < ^ Xfc log fk{l,Qn~k) 



. k=0 



and the following lemma. 



(32) 



Lemma 2 Consider conditional moments of the entity Wn defined at the be- 
ginning of Section 

n-l 

S^:^=Y,M,k, « = 1,2. 

k=0 

Under conditions ^F^, (HI]), and we have 



> X 



Kq 

log a: 



, X — )• oo, 
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and for any fixed e > 0, see (|25p . 



, n — >■ oo. 



Proof of Lemma [5] . For any fixed e > we have 
P (S"^^' > a;) > P ^5*^^^ >x;T< ^mm^ 6li,fe > x"^^ 



> P(St> x^+') - P (T > - P min Oi k < x' 



Notice that according to (|T9|) 

P (T > (loga:)2+'') = o \ t—] , x ^ oo, for any fixed £ > 0. (33) 
\logxJ 

Thus, using ^ and we get 

liminf (log a;- P fs-j!^ > a^H > hminf jloga; • P (5t > x^+')] > Ko/{l + e). 
To obtain a similar estimate from above we write, recall (j25p . 
P f S'i^^ > a;") < P I S'i!^ > x; T < max 6*1 ^ < 

V / \ ^ 0<k<T ' 



+ P (T > Xa) + P \ T < Xc,; max Oi k > x' 

\ 0<k<T ' 

< P(St> x^-"") +P(T > Xa)+P ( max Oi k > x"" 
which together with and ^ yields 

limsup jlogx • P > a;) I < limsup {logx • P (5t > a;-^"^) } 

< i^o/(l-e). 

Finally, according to 



logn 



, n — > oo. 



□ 



Proof of (|30|) . It follows from (|32l) and monotonicity of (3„ that for any fixed 
ee (0,1) 



P (Z„ > 0) > E 



'T-l 



1 - exp <^ ^ Xfe log fk{l, Qn) \ ; 4^ < n'S^j^\T < Uo, 



. k=0 
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Observe that log(l — x) < ~x and 

which holds due to monotonicity of the second derivative of the generating 
function. Therefore, 

log/(l,s)<-(?i(l-s) + (02/2) (1-5)2 
and it follows from ^ that g 



T-l 



T 



fe=0 



for sufficiently large n. As a result we see that for large n 

P (^n > 0) > E [l - e^^^'^r' ; ^ neg{l)^j, < 



1 - e 



' (T > n„) - P (S-^^^ > T < 



In view of ([33]) and Lemma [5] it remains to observe that for any positive A 



E 



1 - e-^^^ 



A ^ 0, 



log(l/A)^ 

which again due to Lemma [2] follows from the Tauberian Theorem 4 in flT, Ch. 
XIIL5] applied to the right hand side of 



A^^E 



1 _ e~^^T 



P ( 5^^^ > a; ) e^^^'da;. 



□ 



Proof of ([3T|) . In a similar manner the inequality 



P {Zn > 0) < E 



'T-l 



1 - exp ^ ^ Xfc log/fc(l, Qn-r) 



. k=0 



entails that for all sufficiently large n 

1 - exp < y^Xk log /fc(l, Qn-r) \\T <na\ max Oi^k < 



E 



. fc=0 



< E 

< E 

< E 



'T-l 



1 - exp <^ ^ Xfc log(l - c0i^kn^^) y,T <na] max^0i^k < 

(T-l 

1 - exp ^ Xfc log(l 



. k=0 



1 - e" 
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Thus 

P (Z„ > 0) < E 



P (T > rio 



P I max Ui k > 

. 0<k<nc 



and (ini) follows due to ^ and 



□ 



5 The subcritical case with an IID environment 

We continue studying BGW-processes in IID environment but now we assume 

E[C]<0, Var[C]G(0,^) (34) 

instead of (|17p . Results rely upon a theorem from [16 giving the asymptotics 
for P(Wt > a;) as X — t' oo. It requires the important technical assumption of 
existence of a constant k such that 

E [e'^'^] = E [nl] = 1, < K < oo. (35) 

If, in addition, for some (5 > 

< E [^2^+''] < oo, E [e^] < oo, (36) 

and 

^>i, E[iein <oo}u{/«<i, E[iM2-M?r + ie2-e?r] <oo}, (a?) 

then, according to 16J, there exists a constant G (0,oo) such that 

P (Wt > x) - C^a;-", x ^ oo. (38) 
Theorem 3 // conditions ([34]), jaS]), ^ hold, then 

P(X„ + Z„>0)-P(Z„>0)- if, -9,(71), n^oo, (39) 
for some positive constant K^, given by (|4ip below, where 



n if K < I, 

q^in) = { n^Mogn, ifn^l, 



(40) 



1 , if K > 1. 

Proof Starting from the first equality in (|32p we obtain for sufhciently large n 



P (Z„ > 0) > E 



l-'[[Ql';T<Nlogn 



k=0 



> E [l 



^Wt log Q, 



] - P(T > iVlogn), 



13 



where a positive constant N will be specified later. On the other hand, we have 
a similar upper bound 

P{Zn>0) < e\i~Q^J^;T <Nlogn\+P{T> Nlogn) 
< E [1 - e^^^°s'5"-"""="] + P {T > Nlogn) . 

Now observe that due to ^ and the same Tauberian theorem used in 
the critical case yields 



with 



^(l-K)c.(4)^ if K<i, 

if K — 1, 



^^^^ 



[ ^J^°^P{WT>x)dx, if K>1. 

It remains to use the fact that (see [4], [5], [6]) under (p4|) and ([35 

P (T > n) = o (A") for some constant A e (0, 1), 
if we choose N = 2/log(A"^). 



(41) 



(42) 

□ 



6 The critical case with a Markovian environ- 
ment 

As compared to the IID case. Markovian environments require extra conditions 
on the underlying Markov chain. In particular, we assume (jl6p and that for 
some p > 

P{t>x)^o( \^ ) , x^oo. (43) 
V X log x/ 

This implies that a :— E[t] < oo and due to conditions E [C*] = and 
E [C] = become equivalent. 

Moreover, under condition (^5)) the sequence of regeneration times (fT^ sat- 
isfies 

P ( \k-^Tk -a\>e) =0 {{logk)-^-P) , k~^oo, (44) 
for an arbitrarily small e > 0, cf. [Tl] . 

Theorem 4 Assume ([H]), (gSl) and 

E[C] = 0, Var[C*] e (0,oo), (45) 

E ^2 O^i)^^ (1 + niax(0,log^*)) < cx), E (mi) < oo. 
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Let, further, for some positive a 

p{ei < i/x)^o{{\ogxr'-") 

P(6»* > x) ^o{{\ogx)-^-") 



> x9l) = o ((log a;) 



-3-a\ 



X — > oo, 
X oo, 
X — > oo. 



Then there exists a constant K* > such that 



P (X„ + Z„ > 0) ~ P (Z„ > 0) 



K* 
logn 



n oo. 



(46) 
(47) 
(48) 



(49) 



Proof The statement is derived in two steps: first 



p {x; + z; > 0) - p (z; > o) 



K* 
logr 



, r — cx) 



(50) 



and then (|49|1 . 

Fix S e (0, 1/4) and write 

P (Z; > 0) = P (Z; > 0; t^s < r^^, \Tr. - ra\ < r^'^) 

+0 (P {t^s > r^^)) + O (P {\Tr ~ ra\ > r6)) . 

Here the last two terms are treated with the help of P (r^a > r^^) < ar~^ and 
(03]), while the main term is analyzed by means of ideas from the proof of 
Theorem [TJ Letting be the number of type 2* daughters produced by 
particles of type 1* and putting T* ~ min{r : Z* = 0}, we deduce from 

P (Z; > 0; T^s < r^\ \Tr - ra\ < r^-^) 

r-l 



E 



1 - n Qr:-r,;rrs < r^' , |r, - ra\ < r^~ 



fc=0 



a lower bound 



P (Z; > 0; T^s < r^\ \Tr -ra\<r 

T*-l 

> E 1 1 - n Q'ri-r, ; T* < /, T,. < r'', \Tr - ra\ < r' 



-<5 



> E 



E 



fe=0 
T'-l 



1- n fkH'^^Q2ra);T* <r\T,s <r^\\Tr-ra\<r^- 



k=Q 

T'-l 



1- n /^(i,Q2™) 



fe=0 



0(P(T*>/))+o — T 



1 



log'+P r 



Hence, applying arguments used to demonstrate (j30p in Lemma [T] one can show 
that 

liminf {logr • P {Z* > 0)} > K* . 
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By similar arguments one can show that 

hmsup {logr ■P{Z;>0)}<K* 



proving ()50|) in view of 



P(X; >0) = ( ^ ) ,r^oo. 



/r 



To demonstrate that follows from (|50l) we set 
N(n) :— max{fc : < n} . 

Clearly, 

P + l + Zkn) + 1 > O) < P (Xn + ^„ > 0) < P [X*^^,,) + > o) , 

(51) 

and for any e E (0, 1) 

P {^Nin) + Z*^in) > O) = P [X*^^,,^ + > 0; N{n) > a-'nil - e)) 

+ 0(P {N{n) < a"in(l -£))) . 

For the first we have 

P {x*Nin) + Z*Nin) > 0; N{n) > a-in(l - e)) < P (x:_,„(i_,) + Z:_,„(i_,) > o) 
On the other hand, again by (|33|) as n — >■ oo 

1 



logi+^n 



P (7V(n) < a-in(l - e)) = P ( 5„-i„(i_,) > n) = 
Thus, 

lim sup {log n ■ P (X„ + Z„ > 0)} 

< limsup {logn . P + Z:-.„(i_,) > O)} < K* . 

A similar estimate from below follows from 

P (-'^Af(Ti) + l + ■^A'(n) + 1 > O) 

> P {K-^mi+e) + K-^mi+e) > O; Nin) + l < a-'nil + s)) 

= P ( K-^n{l+e) + K-^n{l+e) > o) + O ( _ 
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□ 



7 Subcritical processes with a Markovian envi- 
ronment 



Assume now that the two- type subcritical process (X„, Z„) evolves in a station- 
ary Markovian random environment as defined in Section 13.21 Here, similarly 
to Section [S] the auxiliary branching process {X*, Z*) in IID environment with 
probability generating functions (|13p and ([T?l) plays an important role. 

Theorem 5 Assume that 

E[C]<0, Var[C*] e (0,oo), (52) 

and conditions (j35p . (I36p . (|37p are valid for the corresponding random variables 
related to the embedded process (X*, Z*) with the key constant k replaced by k*. 
Suppose, in addition, that 

P (r > x) = o (^2;-i-i«in(«M)^ , x^oo. (53) 

Then there exists a constant K* = K*(k*) > such that, as n — >■ oo, see (|40p . 

P {Xn + Zn > 0) - K*q,, (n) . (54) 



Proof Our main arguments here are similar to that used in the proof of The- 
orem [H Fix £ G (0, 1) and a sufficiently large N and with a = E [r] write 

P(Z;>0) = P{Z; > 0:B{r,e)) 

+0 (P (rNiogr > r"' log^ r) + P (|t^ - ra\ > re)) . 

where 

B (r, e) := ^TNiogr < r'^' log^ r, jr^ - ra\ < erj . 

Clearly, 

P (tn log r > r"* log^ r) < = o (r-^* ) . (55) 

\ / log r V / 

Further, if k* < 1 then, according to [13] we have under condition ([55)) . 

P( -ra| > er) = ofr-"') . (56) 
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Thus, 

F{Z; > 0) > P (Z; > 0;B{r,e)) 



= E 

> E 

> E 

= E 



k=0 
T*-l 

1- n Qr:-r,;T* <N\ogr;B{r,e) 



k=Q 
T'-l 



1- n Q'ra+2er;T* <Nlogr;Bir,e) 



k=0 



I _ logQ,.a + 2£,- 



P(T* > Nlogr) ~T> {B{r,e)) 



where B{r,e) is the event complementary to B{r,e). Due to (fT5)) and we 
have 

P {X; > 0) = P (T* > r) = o {A'') for some ^ < 1. (57) 
It follows that in view of dTSl). (1551) and (1551) one can find iV such that 



P(r* >7Vlogr) + P(S(r,e)) =o 

while using ([2]) and 

P {W^ >y)^ C*y-''' ,C* = C* («*) G (0, oo) , 
one can show, arguing as in Theorem [3l that for k* < 1 



liminf lim r"^ E 



= lim inf lim 



eiO r^oa - log Qra+2Er 



/•oo 

/ P(W^^ > x)e^'°sQ™+2..^2, 



liminf r(l - K*)C , 
eia \TO2a(l + 2e) 



giving for k* < 1 

liminf r'^'P (Z,* > 0) > K^, = T (1 - k*) C 

r— >-oo 

A similar upper bound in view of and §7} yields 

lim r'^'p {x; + z; >0)^K^. 

r—^oo 

If K* > 1 then condition ([5^ entails 

P ( \Tr — ra\ > er) ~ o {r~^^ 



/ 2 Y 
\m2a 
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and, as before, this implies 

lini {q^, (r))-' p {x; + z; > 0) ^ K,, 

r^oo 

where 

= ^^ if .c* = i, 

m,2a \ P {Wf > x) dx, if K* > 1. 
We proceed by recalling ([STj). Now for any e G (0, 1) we have 

P {x*Nin) + Z*^in) > O) = P (X^(„) + Z*^^^) > 0; Nin) > a^'nil - e)) 

+ 0(P ( N{n) < a-^n{l-e))) , 

and as n — > DO 

P (xjV(„) + Z]V(„) > 0; N[n) > a-'n{l - e)) 

It follows from [M] and our conditions that 

P {N{n) < a~^n{l - e)) - P {Sa-in{i-e) > n) = o (n^ ^-^Mk' , n oo. 
Thus, 

hmsup («))"^ P (X„ + Zn > 0) 

< limlimsup {q^,'{n)y^ P {X^(„) + Zj^(„) > Oj 

SO that 

limsup (q.'in))-^ P (X„ + Z„ > 0) < K* , K* = a™^^'"*);^,. 

n— >C30 

The corresponding lower bound is obtained similarly. 

□ 
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